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Abstract 



I^H \ We study the role of collective surface excitations in the electron relaxation 



in small metal particles. We show that the dynamically screened electron- 
electron interaction in a nanoparticle contains a size-dependent correction 



> 
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o. 

^^ , induced by the surface. This leads to new channels of quasiparticle scattering 

accompanied by the emission of surface collective excitations. We calculate 

--^. \ the energy and temperature dependence of the corresponding rates, which 

depend strongly on the nanoparticle size. We show that the surface-plasmon- 

mediated scattering rate of a conduction electron increases with energy, in 

q \ contrast to that mediated by a bulk plasmon. In noble-metal particles, we find 



that the dipole collective excitations (surface plasmons) mediate a resonant 



scattering of ci-holes to the conduction band. We study the role of the latter 

effect in the ultrafast optical dynamics of small nanoparticles and show that, 

with decreasing nanoparticle size, it leads to a drastic change in the differential 

absorption lineshape and a strong frequency dependence of the relaxation near 

the surface plasmon resonance. The experimental implications of our results 

in ultrafast pump-probe spectroscopy are also discussed. 
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I. INTRODUCTION 

The properties of small metal particles in the intermediate regime between bulk-like 
and molecular behavior have been a subject of great interest recently.EfQ Even though the 
electronic and optical properties of nanoparticles have been extensively studied, the effect of 
confinement on electron dynamics is much less understood. Examples of outstanding issues 
include the role of electron-electron interactions in the process of cluster fragmentation, 
the role of surface lattice modes in providing additional channels for intra-molecular en- 
ergy relaxation, the influence of the electron and nuclear motion on the superparamagnetic 
properties of clusters, and the effect of confinement on the nonlinear optical properties and 
transient response under ultrafast excitation.ErBel These and other dynamical phenomena can 
be studied with femtosecond nonlinear optical spectroscopy, which allows one to probe the 
time evolution of the excited states with a resolution shorter than the energy relaxation or 
dephasing times. 

Surface collective excitations play an important role in the absorption of light by metal 
nanoparticles. In large particles with sizes comparable to the wave-length of light A (but 
smaller than the bulk mean free path), the lineshape of the surface plasmon (SP) resonance 
is determined by the electromagnetic effects.^ In small nanoparticles with radii i? <C A, 
the absorption spectrum is governed by quantum confinement effects. For example, the 
momentum non-conservation due to the confining potential leads to the Landau damping 
of the SP and to a resonance linewidth inversely proportional to the nanoparticle size.ErB 
Confinement changes also non-linear optical properties of nanoparticles: a size-dependent 
enhancement of the third order susceptibilities, caused by the elastic surface scattering of 
single-particle excitations, has been reported.Qu 

Extensive experimental studies of the electron relaxation in nanoparticles have recently 
been performed using ultrafast pump-probe spectroscopyETEl Unlike in semiconductors, the 
dephasing processes in metals are very fast, and nonequilibrium populations of optically 
excited electrons and holes are formed within several femtoseconds. These thermalize into 



the hot Fermi-Dirac distribution within several hundreds of femtoseconds, mainly due to 
e-e and h-h scatteringJlZTEj Since the electron heat capacity is much smaller than that of 
the lattice, a high electron temperature can be reached during less than 1 ps time scales, 
i.e., before any significant energy transfer to the phonon bath occurs. During this stage, the 
SP resonance was observed to undergo a time-dependent spectral broadeningJiirc^l Subse- 
quently, the electron and phonon baths equilibrate through the electron-phonon interactions 
over time intervals of a few picoseconds. During this incoherent stage, the hot electron dis- 
tribution can be characterized by a time-dependent temperature. Correlation effects play an 
important role in the latter regime. For example, in order to explain the differential absorp- 
tion lineshape, it is essential to take into account the e-e scattering of the optically-excited 
carriers near the Fermi surface.EJ Furthermore, despite the similarities to the bulk-like be- 
havior, observed, e.g., in metal films, certain aspects of the optical dynamics in nanoparticles 
are significantly different .ESLilO For example, experimental studies of small Cu nanoparti- 
cles revealed that the relaxation times of the the pump-probe signal depend strongly on 
frequency: the relaxation was considerably slower at the SP resonance. £3>$3 This and other 
observations suggest that collective surface excitations play an important role in the electron 
dynamics in small metal particles. 

Let us recall the basic facts regarding the linear absorption by metal nanoparticles em- 
bedded in a medium with dielectric constant e m . We will focus primarily on noble metal 
particles containing several hundreds of atoms; in this case, the confinement affects the ex- 
tended electronic states even though the bulk lattice structure has been established. When 
the particles radii are small, R <C A, so that only dipole surface modes can be optically 
excited and non-local effects can be neglected, the optical properties of this system are 
determined by the dielectric function!^ 

e co i{uJ) = e m + 3pe m —-——- — , (1) 

e(w) + 2e m 

where e{u>) = e'(uj) + ie"(u) is the dielectric function of a metal particle and p <C 1 is 
the volume fraction occupied by nanoparticles in the colloid. Since the <i-electrons play an 



important role in the optical properties of noble metals, the dielectric function e(u>) includes 
also the interband contribution q(cj). For p <C 1, the absorption coefficient of such a system 
is proportional to that of a single particle and is given byS 

a(u) = -9pe^-lm^- (2) 

c e s {u) 

where 

e s (u) = e d (uj) -u 2 p /uj(uj + i^ s ) + 2e m , (3) 

plays the role of an effective dielectric function of a particle in the medium. Its zero, 
e' s (uj s ) = 0, determines the frequency of the SP, uj s . In Eq. (||), u p is the bulk plasmon 
frequency of the conduction electrons, and the width 7 S characterizes the SP damping. The 
semiclassical result Eqs. (|2]) and (Q) applies to nanoparticles with radii R ^> q~ x , where q TF 
is the Thomas-Fermi screening wave-vector (q~ F ~ 1 A in noble metals). In this case, the 
electron density deviates from its classical shape only within a surface layer occupying a small 
fraction of the total volume£ll Quantum mechanical corrections, arising from the discrete 
energy spectrum, lead to a width 7 S ~ v F /R, where v F = k F /m is the Fermi velocity.El'Q 
Even though ^ s /uj s ~ (q TF R)~ 1 <C 1, this damping mechanism dominates over others, e.g., 
due to phonons, for sizes R ^ 10 nm. In small clusters, containing several dozens of atoms, 
the semiclassical approximation breaks down and density functional or ab initio methods 
should be used.ElQ 

It should be noted that, in contrast to surface collective excitations, the e-e scattering 
is not sensitive to the nanoparticle size as long as the condition q TF R 3> 1 holds.E^ Indeed, 
for such sizes, the static screening is essentially bulk-like. At the same time, the energy 
dependence of the bulk e-e scattering rate,c3 ^ e oc (E—Ep) 2 , with Ep being the Fermi energy, 
comes from the phase-space restriction due to the momentum conservation, and involves the 
exchange of typical momenta q ~ q TF . If the size-induced momentum uncertainty 5q ~ R~ x 
is much smaller than q TF , the e-e scattering rate in a nanoparticle is not significantly affected 
by the confinement." 



In this paper we address the role of collective surface excitations in the electron re- 
laxation in small metal particles. We show that the dynamically screened e-e interaction 
contains a correction originating from the surface collective modes excited by an electron 
in nanoparticle. This opens up new quasiparticle scattering channels mediated by surface 
collective modes. We derive the corresponding scattering rates, which depend strongly on 
the nanoparticle size. The scattering rate of a conduction electron increases with energy, in 
contrast to the bulk-plasmon mediated scattering. In noble metal particles, we study the 
SP-mediated scattering of a <i-hole into the conduction band. The scattering rate of this 
process depends strongly on temperature, and exhibits a peak as a function of energy due 
to the restricted phase space available for interband scattering. We show that this effect 
manifests itself in the ultrafast nonlinear optical dynamics of nanometer-sized particles. In 
particular, our self-consistent calculations show that, near the SP resonance, the differential 
absorption lineshape undergoes a dramatic transformation as the particle size decreases. We 
also find that the relaxation times of the pump-probe signal depend strongly on the probe 
frequency, in agreement with recent experiments. 

The paper is organized as follows. In Section [IT] we derive the dynamically screened 
Coulomb potential in a nanoparticle. In Section [TTIJ we calculate the SP-mediated quasipar- 



ticle scattering rates of the conduction electrons and the c?-band holes. In Section [TV] we 
incorporate these effects in the calculation of the absorption spectrum and study their role 
in the size and frequency dependence of the time-resolved pump-probe signal. 

II. ELECTRON ELECTRON INTERACTIONS IN METAL NANOPARTICLES 

In this section, we study the effect of the surface collective excitations on the e-e inter- 
actions in a spherical metal particle. To find the dynamically screened Coulomb potential, 
we generalize the method previously developed for calculations of local field corrections to 
the optical fields.^ The potential U(u; r, r') at point r arising from an electron at point r' 
is determined by the equation^ 



U(u; r, r') = u(r - r') + / dr l dr 2 u(r - ri)n(w; ri, r 2 )E/(w; r 2 , r'), (4) 

where w(r — r') = e 2 |r — r'| _1 is the unscreened Coulomb potential and n(c<j; ri, r 2 ) is the 
polarization operator. There are three contributions to II, arising from the polarization of 
the conduction electrons, the <i-electrons, and the medium surrounding the nanoparticles: 
II = II C + ILj + IT m . It is useful to rewrite Eq. (f|) in the "classical" form 

V- (E + 4ttP) = 47re 2 <5(r-r'), (5) 

where E(u>; r, r') = — V?7(u>; r, r') is the screened Coulomb field and P = P c + P^ + P m is 
the electric polarization vector, related to the potential U as 

VP(cj; r, r') = -e 2 f dr^u; r, r x )U(uj; r x , r'). (6) 

In the random phase approximation, the intraband polarization operator is given by 

n c (c,;r,rO =E i ( ^ c ~f ^In CW^WC^)^^), (7) 

aa' L a ~ &a> + U + t[) 

where E c a and ip%, are the single-electron eigenenergies and eigenfunctions in the nanopar- 
ticle, and f(E) is the Fermi-Dirac distribution (we set h — 1). Since we are interested in 
frequencies much larger than the single-particle level spacing, II c (u;) can be expanded in 
terms of l/u>. For the real part, 11^ (a;), we obtain in the leading orderEl 

IT>; r, n) = ^VK(r)V5(r - n)], (8) 

where n c (r) is the conduction electron density. In the following we assume, for simplicity, 
a step density profile, n c (r) = n c 9(R — r), where n c is the average density. The leading 
contribution to the imaginary part, U"(u>), is proportional to u~ 3 , so that 11" (a;) <C U' c (uj). 
By using Eqs. (|8]) and (§), one obtains a familiar expression for P c at high frequencies, 

P c (u; r, r') = £^W(u;; r, r') = 6(R - r) Xc (u;)E(u;; r, r'), (9) 

where Xc(w) = —e 2 n c /muj 2 is the conduction electron susceptibility. Note that, for a step 
density profile, P c vanishes outside the particle. The <i-band and dielectric medium contri- 
butions to P are also given by similar relations, 



P d (uj; r, r') = 6{R - r) Xd (u;)E(u;; r, r'), (10) 

P m (u; r, r') = 9(r - i?)x m E(u;; r, *'), (11) 

where Xi — ( e i ~~ 1)/47T, i = d,m are the corresponding susceptibilities and the step functions 
account for the boundary conditions.B3 Using Eqs. ©-(11]), one can write a closed equation 
for U(u>; r, r'). Using Eq. @, the second term of Eq. (Q) can be presented as — e~ 2 J driu(r — 
ri) V • P(o>; ri, r'). Substituting the above expressions for P, we then obtain after integrating 
by parts 

-,2 



e(u)U(u; r, r') = -^r+ f ^Vi^- • Vi [9(R - r)x(w) + 0(r - fl)xm] l/(u>; r 1? r') 

|r — r | J |r — ri| 

+i f dv l dv 2 —^ — r n' c '(u;; ri, r 2 )J7(w; r 2 , r'), (12) 

J r — rd 



; 2 

|r — ri 
with 



e(w) = 1 + 4ttxM = e d (w) - wj/w 2 , (13) 

a; 2 = 47re 2 n c /m being the plasmon frequency in the conduction band. The last term in 
the rhs of Eq. ([T2]), proportional to H"(u), can be regarded as a small correction. To solve 
Eq. (|T2"D, we first eliminate the angular dependence by expanding U(u;r,r') in spherical 
harmonics, Yi,m{t), with coefficients ULM(u;r,r'). Using the corresponding expansion of 
|r — r'| _1 with coefficients QLhiij"^') = 2L+i r ~ L ~ lr ' L (for r > r'), we get the following 
equation for Ulm(w, r,r') : 

L + 1 / t \ ^ 
e{u)U LM {ix)] r, r') = Qlm(t, r') + 4tt [x(u;) - x™] f -= J ^lm(w; #, r') 

+^ 2 H dr 1 dr 2 r 2 l rlQ LM {r,r l )Ill ML , MI (uj;r 1 ,r2)U L 'M'{uJ;r2,r'), (14) 

where 

n LM,L'Af'( w ; r i,r 2 ) = / dr 1 dr 2 Yl M (r 1 )U"(uj;r 1 ,r2)Y LI M'(i2), (15) 

are the coefficients of the multipole expansion of II" (u;; ri, r 2 ). For 11" = 0, the solution of 
Eq. ( |14"D can be presented in the form 



47Te 2 rpL^rL 

U LM {u- r, r') = a(u)e 2 Q LM (r, r') + b(u) ^ - ^ R2L+l , (16) 

with frequency-dependent coefficients a and b. Since Tl"(u>) <C n^(a>) for relevant frequen- 
cies, the solution of Eq. (|T4] ) in the presence of the last term can be written in the same 
form as Eq. (|T6"D, but with modified a(u) and b{uS). Substituting Eq. ( |I6"D into Eq. ([14]), we 
obtain after lengthy algebra in the lowest order in II" 

a(u;)=e- 1 (uj) 1 b(oo) = e L 1 (uj)-e- 1 (uj), (17) 

where 

Cl( ^ = 2L + T e{uj) + 2LTl em + ^ l( ^' (18) 

is the effective dielectric function, whose zero, e' L {uJi) = 0, determines the frequency of the 
collective surface excitation with angular momentum L,0 



Lu, 



Ul 



(19) 



L Le' d (u L ) + (L + l)e m 
In Eq. (|T8|), e" L (u;) characterizes the damping of the L-pole collective mode by single-particle 
excitations, and is given by 

4-7T-V _ 

<M = {2L - 1)jR2L+1 E \M L J\ 2 [f{El) - f(E°,)]5(K - K, + u), (20) 

where M^f are the matrix elements of t l Ylm (r). Due to the momentum nonconservation 
in a nanoparticle, the matrix elements are finite, which leads to the size-dependent width 
of the L-pole mode:BEJ 

1l = 2J ^- 2 eM- (21) 

For oj ~ u L , one can show that the width, ^ L ~ Vp/R, is independent of to. Note that, in 
noble metal particles, there is an additional (i-electron contribution to the imaginary part 
of 6l(^) at frequencies above the onset A of the interband transitions. 

Putting everything together, we arrive at the following expression for the dynamically- 
screened interaction potential in a nanoparticle: 

8 



J\ „2 A„ i /„„. A L 



^-■^-^^S^T^^j r^)n^), (22) 

with 1~[}{uj) = ej}{ui) — e _1 (c<j). Equation fl22|), which is the main result of this section, 
represents a generalization of the plasmon pole approximation to spherical particles. The 
two terms in the rhs describe two distinct contributions. The first comes from the usual 
bulk-like screening of the Coulomb potential. The second contribution describes a new 
effective e-e interaction induced by the surface: the potential of an electron inside the 
nanoparticle excites high-frequency surface collective modes, which in turn act as image 
charges that interact with the second electron. It should be emphasized that, unlike in the 
case of the optical fields, the surface-induced dynamical screening of the Coulomb potential 
is size-dependent. 

Note that the excitation energies of the surface collective modes are lower than the bulk 
plasmon energy, also given by Eq. ([19]) but with e m = 0. This opens up new channels of 
quasiparticle scattering, considered in the next section. 

III. QUASIPARTICLE SCATTERING VIA SURFACE COLLECTIVE MODES 

In this section we calculate the rates of quasiparticle scattering accompanied by the 
emission of surface collective modes. We start with the scattering of an electron in the 
conduction band. In the first order in the surface-induced potential, given by the second 
term in the rhs of Eq. (|22|), the corresponding scattering rate can be obtained from the 
Matsubara self-energyc3 

E - M = -fe £? wi^ l M G ' A -' + *">• (23) 

where G c a = (iou — E^.) -1 is the non-interacting Green function of the conduction electron. 
Here the matrix elements M^£f are calculated with the one-electron wave functions ^,(r) = 
Rni(r)Yi m (r) . Since \a) and \a') are the initial and final states of the scattered electron, the 
main contribution to the Lth term of the angular momentum sum in Eq. (^) will come from 
electron states with energy difference E a — E a / ~ u>l- Therefore, M^I can be expanded 



in terms of the small parameter Eq/\E^ — E c a ,\ ~ Eq/ujl, where Eq = (2mR 2 )~ 1 is the 
characteristic confinement energy. The leading term can be obtained by using the following 
procedure. Sc3 We present M^f as 

MSf = <e,a|r^ M (r)|e, a <> = MlgJ^M!^ , (24) 

where H = H Q + V(r) is the Hamiltonian of an electron in a nanoparticle with confining 
potential V(r) = Vq9(t — R). Since [H,r L Yi,M{?j\ = — ^V[r L Yi,M(r)] • V, the numerator in 
Eq. (|24l) contains a term proportional to the gradient of the confining potential, which peaks 
sharply at the surface. The corresponding contribution to the matrix element describes the 
surface scattering of an electron making the L-pole transition between the states \c,a) 
and \c, a'), and gives the dominant term of the expansion. Thus, in the leading order in 
\E^ — £J°/| -1 , we obtain 

m lm (c,a|V[r^ LM (f)]-W(r)|c,a') LR L ^ ^ LM 

with (pfJf Vm i = J dYY* m (Y)Y LM {i)Yii m i{v). Note that, for L = 1, Eq. (ESI) becomes exact. For 



electron energies close to the Fermi level, E^ t ~ Ep, the radial quantum numbers are large, 
and the product V R n i(R)R n 'i'(R) can be evaluated by using semiclassical wave-functions. 
In the limit V — »■ oo, this product is given byB 2JE^ ll E^, l ,/R 3 , where E c nl = n 2 (n + l/2) 2 E 
is the electron eigenenergy for large n. Substituting this expression into Eq. (|25|) and then 
into Eq. (0), we obtain 

£«M = ~ E E E £4 t^t; 7#%ti S G «' (ay + !w) ' (26) 

with 

^ = E Kiwi 2 = (2L + ^f + 1} /* dxPi{x)P L (x)M*), (27) 

where Pj(a;) are Legendre polynomials; we used properties of the spherical harmonics in the 
derivation of Eq. (ETj). For E^ ~ £7?, the typical angular momenta are large, I ~ k F R ~^> 1, 
and one can use the large-/ asymptotics of Pi; for the low multipoles of interest, L <C I, 

10 



the integral in Eq. ( p7|) can be approximated by qF+i^ 11 '' After performing the Matsubara 
summation, we obtain for the imaginary part of the self-energy that determines the electron 
scattering rate 

where N(E) is the Bose distribution and gi(E) is the density of states of a conduction 
electron with angular momentum I, 



gi (E) = 2^6(E° nl -E)^^ 2 -^, (29) 

where we replaced the sum over n by an integral (the factor of 2 accounts for spin) . 

Each term in the sum in the rhs of Eq. fl2"8|) represents a channel of electron scattering 
mediated by a collective surface mode with angular momentum L. For low L, the difference 
between the energies of modes with successive values of L is larger than their widths, so 
that the different channels are well separated. Note that since all ujl are smaller than the 
frequency of the (undamped) bulk plasmon, one can replace €l{u) by €l(u) in the integrand 
of Eq. (]28p for frequencies u> ~ u^. 

Consider now the L = 1 term in Eq. fl28|) , which describes the SP-mediated scattering 
channel. The main contribution to the integral comes from the SP pole in e[ {oj) = 3e~ 1 (c<j), 
where e a (u) is the same as in Eq. (||). To estimate the scattering rate, we approximate 
Ime7 1 (co') by a Lorentzian, 

%uj 2 p /uj 3 + e" d {u) u 2 s u sl 



Ime c (u>) 



[e'H + 2e m ] 2 + [7,^/^3 + e ^ u)] 2 ^j + 2em (w a _ w a)2 + ^ > 

(30) 



where ui s = U\ = U) p /Je d (u> s ) + 2e m and 7 = 7 S + iO s e d {u s ) are the SP frequency and width, 
respectively. For typical widths 7 <C u s , the integral in Eq. fl2"g| ) can be easily evaluated, 
yielding 



24e 2 co>..E , n E c n J2m(uj - u s ) 
Im£° (u) = - f eu} **° "V ^ J^h _ f , _ )i (31) 

al ; e^w.) + 2e m (w - El - u s y L M S/M l ; 

11 



Finally, using the relation e 2 fcf[e^(o; s ) + 2e m ] 1 = 3tiuj 2 /8Ef, the SP-mediated scattering 
rate, "/l(E°) = — ImE°(E^), takes the form 

lliE) = 9^ JL (^) V2 [1 - f(E - *.)]. (32) 

Recalling that E = (2mR 2 )~ 1 , we see that the scattering rate of a conduction electron 
is size-dependent: 7f ^ R~ A . At E = E F + cu s , the scattering rate jumps to the value 
97r(l + uj s /Ef)Eq/lu s , and then increases with energy as E 3 ^ 2 (for uj s <C -Ef). This should 
be contrasted with the usual (bulk) plasmon-mediated scattering, originating from the first 
term in Eq. (p2|) , with the rate decreasing as E~ l l 2 above the onset£5l To estimate the size 
at which 7| becomes important, we should compare it with the Fermi liquid e-e scattering 
ratejHa "f e (E) = \^ ~ E J ■ For energies E ~ E F + lu s , the two rates become comparable 
for 



16fc F 



(KK) ^ l2 EL( l + ELfl±^\'\ (33) 

u s \ oj s J \nq TF J 

In the case of a Cu nanoparticle with lo s ~ 2.2 eV, we obtain k F R ~ 8, which corresponds to 
the radius R ~ 3 nm. At the same time, in this energy range, the width j* exceeds the mean 
level spacing 5, so that the energy spectrum is still continuous. The strong size dependence 
of 7| indicates that, although 7J increases with energy slower than j e , the SP-mediated 
scattering should dominate for nanometer-sized particles. Note that the size and energy 
dependences of scattering in different channels are similar. Therefore, the total scattering 
rate as a function of energy will represent a series of steps at the collective excitation energies 
E = oj l < u p on top of a smooth energy increase. We expect that this effect could be 
observed experimentally in time-resolved two-photon photoemission measurements of size- 
selected cluster beams." 

We now turn to the interband processes in noble metal particles and consider the scatter- 
ing of a d-hole into the conduction band. From now on we restrict ourselves to the scattering 
via the dipole channel, mediated by the SP. The corresponding surface-induced potential, 
given by the L = 1 term in Eq. (p2|), has the form 

12 



^p 2 r ■ r' 1 

U -^ = H^1M- < 34 > 

With this potential, the ci-hole self-energy is given by 

s«M = -— £|cW| -£ (iw0 ■ (35) 

where d aa / = (c, a\r\d, a') = (c,a\p\d,a')/im(E^ — E%,) is the interband transition matrix 
element. Since the final state energies in the conduction band are high (in the case of 
interest here, they are close to the Fermi level), the matrix element can be approximated by 
the bulk-like expression (c,a\p\d,a') = S aa i(c\p\d) = 5 aa >fi, the corrections due to surface 
scattering being suppressed by a factor of (k F R)~ l <C 1. After performing the frequency 
summation, we obtain for ImS^ 

i v'r ^ 9g V , jng| zA± /(*%) M 

toE » = - ro . (Jg) . g fa u (,Ei- u ) ' (36) 

with Eff = E c a — E*. We see that the scattering rate of a ci-hole with energy E%, Jhi-^a) = 
Im£^(£^), has a strong R~ 3 dependence on the nanoparticle size, which is, however, different 
from that of the intraband scattering, Eq. fl32|) . 

The important difference between the interband and the intraband SP-mediated scat- 
tering rates lies in their energy dependence. Since the surface-induced potential, Eq. (fill), 
allows for only vertical (dipole) interband single-particle excitations, the phase space for the 
scattering of a d-ho\e with energy E^ is restricted to a single final state in the conduction 
band with energy E c a . As a result, the <i-hole scattering rate, , jf l (E^), exhibits a peak as 
the difference between the energies of final and initial states, E^f = E c a — E^, approaches 
the SP frequency u s [see Eq. fl36|)1. In contrast, the energy dependence of 7! is smooth due 
the larger phase space available for scattering in the conduction band. This leads to the 
additional integral over final state energies in Eq. (|28|) , which smears out the SP resonant 
enhancement of the intraband scattering. 

As we show in the next section, the fact that the scattering rate of a <i-hole is dominated 
by the SP resonance, affects strongly the nonlinear optical dynamics in small nanoparticles. 

13 



This is the case, in particular, when the SP frequency, u> a , is close to the onset of interband 
transitions, A, as, e.g., in Cu and Au nanoparticles.uOoEj Consider an e-h pair with 
excitation energy uj close to A. As we discussed, the ci-hole can scatter into the conduction 
band by emitting a SP. According to Eq. (0), for uj ~ u> s , this process will be resonantly 
enhanced. At the same time, the electron can scatter in the conduction band via the usual 
two-quasiparticle process. For uj ~ A, the electron energy is close to Ep, and its scattering 
rate is estimated asc3 7 e ~ 10 -2 eV. Using the bulk value of fi, 2jj 2 /m ~ 1 eV near the 
L-point,c2l we find that 7^ exceeds 7 e for R < 2.5 nm. In fact, one would expect that, in 
nanoparticles, \x is larger than in the bulk due to the localization of the conduction electron 
wave-functions.u 

IV. SURFACE PLASMON NONLINEAR OPTICAL DYNAMICS 

In this section, we study the effect of the SP-mediated interband scattering on the non- 
linear optical dynamics in noble metal nanoparticles. When the hot electron distribution 
has already thermalized and the electron gas is cooling to the lattice, the transient response 
of a nanoparticle can be described by the time-dependent absorption coefficient a(u,t), 
given by Eq. fl2|) with time-dependent temperature. E3 In noble-metal particles, the temper- 
ature dependence of a originates from two different sources. First is the phonon-induced 
correction to 7 S , which is proportional to the lattice temperature Ti(t). As mentioned in the 
Introduction, for small nanoparticles this effect is relatively weak. Second, near the onset of 
the interband transitions, A, the absorption coefficient depends on the electron temperature 
T(t) via the interband dielectric function 6,1(00) [see Eqs. (0) and @]. In fact, in Cu or Au 
nanoparticles, u s can be tuned close to A, so the SP damping by interband e-h excitations 
leads to an additional broadening of the absorption peakJil In this case, the temperature 
dependence of e<j(u>) dominates the pump-probe dynamics. Below we show that, near the 
SP resonance, both the temperature and frequency dependence of e<j (a>) = 1 + AitXdi^) are 
strongly affected by the SP-mediated interband scattering. 
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For non-interacting electrons, the interband susceptibility, Xd{iw) — Xd(iw) + x<j(— zo;), 
has the standard former 

M**) = -E^|^^E^M)^K + ^), (37) 

where G^(iu') is the Green function of a d-electron. Since the <i-band is fully occupied, the 
only allowed SP-mediated interband scattering is that of the d— hole. We assume here, for 
simplicity, a dispersionless d-band with energy E d . Substituting G^fiu') = [iu' — E d + E F — 
E^(io/)] _1 , with T^iiuS) given by Eq. fl35"D , and performing the frequency summation, we 
obtain 

. eV r dE c g(E c ) f(E c )-l 

Xd{U> m 2 J (E cd ) 2 u-E cd + iY h {u,E c y { ' 

where g{E c ) is the density of states of conduction electrons. Here Jh(u), E c ) = lmS d (E c — u) 
is the scattering rate of a <i-hole with energy E c — u, for which we obtain from Eq. (|3~6|), 



where we neglected N{u>) for frequencies uo ~ uo s ^> k^T . Remarkably, 7^(0;, _E C ) exhibits 
a sharp peak as a function of the frequency of the probe optical field. The reason for this 
is that the scattering rate of a <i-hole with energy E depends explicitly on the difference 
between the final and initial states, E c — E, as discussed in the previous section: therefore, 
for a (i-hole with energy E = E c — u, the dependence on the final state energy, E c , cancels 
out in e s (E c — E) [see Eq. ([36])]. This implies that the optically-excited <i-hole experiences 
a resonant scattering into the conduction band as the probe frequency u approaches the SP 
frequency. It is important to note that 7^(0;, _E C ) is, in fact, proportional to the absorption 
coefficient a(u) [see Eq. (0)]. Therefore, the calculation of the absorption spectrum is a 
self-consistent problem defined by Eqs. (Q), (|3|), (|38|), and (|39|). 

It should be emphasized that the effect of 7^ on e'j(u) increases with temperature. Indeed, 
the Fermi function in the rhs of Eq. (^) implies that 7^ is small unless E c — Ep < ksT. 
Since the main contribution to x'd(w) comes from energies E c — Ep ~ uj — A, the <i-hole 
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scattering becomes efficient for electron temperatures k B T > uj s — A. As a result, near the 
SP resonance, the time evolution of the differential absorption, governed by the temperature 
dependence of a, becomes strongly size-dependent, as we illustrate in the rest of this section. 
In the numerical calculations below, we adopt the parameters of the experiment of Ref. 
|TT| , which was performed oni?~ 2.5 nm Cu nanoparticles with SP frequency, uo s ~ 2.22 eV, 
slightly above the onset of the interband transitions, A ~ 2.18 eV. In order to describe the 
time-evolution of the differential absorption spectra, we first need to determine the time- 
dependence of the electron temperature, T(t), due to the relaxation of the electron gas to 
the lattice. For this, we employ a simple two-temperature model, defined by heat equations 
for T(t) and the lattice temperature T}(£): 

C(T)^ = -G(T-T l ) 1 

BT 
C l -^ = G(T-T l ), (40) 

where C(T) = FT and C\ are the electron and lattice heat capacities, respectively, and G 
is the electron-phonon coupling.!^ The parameter values used here were G = 3.5 x 10 16 
Wm~ 3 K~ x , T = 70 Jm _3 K -2 , and Ci = 3.5 Jm~ 3 K . The values of 7 S and \i were extracted 
from the fit to the linear absorption spectrum, and the initial condition for Eq. ( fJTO|) was 
taken as T = 800 K, the estimated pump-induced hot electron temperature. E-il We then self- 
consistently calculated the time-dependent absorption coefficient a(oo,t), and the differential 
transmission is proportional to a r (u) — a(ui,t), where a r {uj) was calculated at the room 
temperature. 

In Fig. 1 we plot the calculated differential transmission spectra for different nanoparticle 
sizes. Fig. 1(a) shows the spectra at several time delays for R = 5.0 nm; in this case, the 
SP-mediated d-hole scattering has no significant effect. Note that it is necessary to include 
the intraband e-e scattering in order to reproduce the differential transmission lineshape 
observed in the experiment^ For optically excited electron energy elo S e to E F , this can be 
achieved by adding the e-e scattering rateca r ) e {E c ) oc [1 — f(E c )][(E c — E F ) 2 + (irk B T) 2 ] 
to 7^ in Eq. (|38|) . The difference in r ) e {E c ) for E c below and above Ep leads to a lineshape 
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similar to that expected from the combination of red-shift and broadening. 

In Figs. 1(b) and (c) we show the differential transmission spectra with decreasing 
nanoparticle size. For R = 2.5 nm, the apparent red-shift is reduced [see Fig. 2(b)]. This 
change can be explained as follows. Since here uj s ~ A, the SP is damped by the interband 
excitations for u > oo s , so that the absorption peak is asymmetric. The ci-hole scattering 
with the SP enhances the damping; however, since the tu-dependence of 7^ follows that of 
a, this effect is larger above the resonance. On the other hand, the efficiency of scattering 
increases with temperature, as discussed above. Therefore, for short time delays, the in- 
crease in the absorption is relatively larger for 00 > uo s . With decreasing size, the strength of 
this effect increases, leading to an apparent blue-shift [see Fig. 2(c)]. Such a strong change 
in the absorption dynamics originates from the R~ 3 dependence of the <i-hole scattering 
rate; reducing the size by the factor of two results in an enhancement of 7^ by an order of 
magnitude. 

In Fig. 2 we show the time evolution of the differential transmission at several frequencies 
close to uo s . It can be seen that the relaxation is slowest at the SP resonance; this char- 
acterizes the robustness of the collective mode, which determines the peak position, versus 
the single-particle excitations, which determine the resonance width. For larger sizes, at 
which 7^ is small, the change in the differential transmission decay rate with frequency is 
smoother above the resonance [see Fig. 2(a)]. This stems from the asymmetric lineshape 
of the absorption peak, mentioned above: the absorption is larger for uj > u s , so that its 
relative change with temperature is weaker. For smaller nanoparticle size, the decay rates 
become similar above and below u s [see Fig. 2(b)]. This change in the frequency dependence 
is related to the stronger SP damping for u > u s due to the d-ho\e scattering, as discussed 
above. Since this additional damping is reduced with decreasing temperature, the relaxation 
is faster above the resonance, compensating the relatively weaker change in the absorption. 
This rather "nonlinear" relation between the time-evolution of the pump-probe signal and 
that of the temperature, becomes even stronger for smaller sizes [see Fig. 2(c)]. In this 
case, the frequency dependence of the differential transmission decay below and above uu s is 
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reversed. Note, that a frequency dependence consistent with our calculations presented in 



Fig. 2(b) was, in fact, observed in the experiment of Ref. |TT]. At the same time, the changes 
in the linear absorption spectrum are relatively small. 

V. CONCLUSIONS 

To summarize, we have examined theoretically the role of size-dependent correlations in 
the electron relaxation in small metal particles. We identified a new mechanism of quasi- 
particle scattering, mediated by collective surface excitations, which originates from the 
surface-induced dynamical screening of the e-e interactions. The behavior of the corre- 
sponding scattering rates with varying energy and temperature differs substantially from 
that in the bulk metal. In particular, in noble metal particles, the energy dependence of the 
(i-hole scattering rate was found similar to that of the absorption coefficient. This led us 
to a self-consistent scheme for the calculation of the absorption spectrum near the surface 
plasmon resonance. 

An important aspect of the SP-mediated scattering is its strong dependence on size. 
Our estimates show that it becomes comparable to the usual Fermi-liquid scattering in 
nanometer-sized particles. This size regime is, in fact, intermediate between "classical" 
particles with sizes larger than 10 nm, where the bulk-like behavior dominates, and very 
small clusters with only dozens of atoms, where the metallic properties are completely lost. 
Although the static properties of nanometer-sized particles are also size-dependent, the 
deviations from their bulk values do not change the qualitative features of the electron 
dynamics. In contrast, the size-dependent many-body effects, studied here, do affect the 
dynamics in a significant way during time scales comparable to the relaxation times. As 
we have shown, the SP-mediated interband scattering reveals itself in the transient pump- 
probe spectra. In particular, as the nanoparticle size decreases, the calculated time-resolved 
differential absorption develops a characteristic lineshape corresponding to a resonance blue- 
shift. At the same time, near the SP resonance, the scattering leads to a significant change 



in the frequency dependence of the relaxation time of the pump-probe signal, consistent 
with recent experiments. These results indicate the need for a systematic experimental 
studies of the size-dependence of the transient nonlinear optical response, as we approach 
the transition from boundary-constrained nanoparticles to molecular clusters. 
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FIGURES 
FIG. 1. Calculated differential transmission spectra at positive time delays for nanoparticles 

with (a) R = 5 nm, (b) R = 2.5 nm, and (c) R = 1.2 nm. 

FIG. 2. Temporal evolution of the differential transmission at frequencies close the SP reso- 
nance for nanoparticles with (a) R = 5 nm, (b) R = 2.5 nm, and (c) R = 1.2 nm. 
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